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Abstract 

Minimizing the effect of decoherence on a quantum register must be a central part of any strategy 
to realize scalable quantum information processing. Apart from the strength of the coupling to 
the environment, the decoherence rate is determined by the the system level structure and by 
the spectral composition of the noise trace that the environment generates. Here, we discuss a 
relatively simple model that allows us to study these different effects quantitatively in detail. We 
evaluate the effect that the perturbation has on a nuclear magnetic resonance (NMR) system while 
it performs a Grover search algorithm. 
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I. INTRODUCTION 



The advantage that quantum computers have over classical computers hinges on the 
creation and preservation of quantum coherence Any real quantum computer, however, 
interacts with its environment and such interactions result in decoherence which increases 
the probability that the quantum computation may fail [2]. Decoherence is thus one of the 
main obstacles for building practical quantum computers. 

A number of strategies have developed for suppressing decoherence, including quantum 
error corrections [3j], dynamical decoupling (quantum control) [4|, decoherence-free subspaces 
0, holono m ic — computation & the quantum Zeno effect g and spectra, degen- 
eracy systems [8|, [9( . It has been proved that the first four strategies can be unified under 
a general algebraic framework |lOj, and the quantum Zeno effect can be unified with dy- 



namical decoupling [ll|]. Up to the present quantum error corrections [12j|, decoherence-free 



subspaces 



13|, holonomic quantum computation I4J], and the suppression of artificial deco 



herence by dynamical decoupling (bang-bang control) [15| have been experimentally tested 
using nuclear magnetic resonance (NMR). 

In this paper, we concentrate on a different aspect: The decoherence rate is not only 
determined by the strength of the coupling to the environment, but also by its operational 
form. As an example, the decoherence differs qualitatively if the coupling operator commutes 
with the system operator. Another aspect is the spectral composition of the noise: If 
the environment is (almost) static, the interaction is adiabatic. If it has components that 
fluctuate at transition frequencies of the system Hamiltonian, its effect can be particularly 
strong. 

A number of model systems have been discussed to study the interaction of a quantum 
register with a noisy environment. One model is known as the spin bath where the envi- 
ronment consists of a set of two- level systems or spin- 1/2 systems 16) . In another model, 
;he so-called spin-boson model, the environment consists of a set of harmonic oscillators 
NMR can simulate the decoherence effect (or called artificial decoherence) 
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22J or imple- 



through the interactions generated by the spins viewed as spin bath 
ment error models by radio- frequency and gradient pulses for demonstrating quantum error 
corrections 

Here, we use a semiclassical model, where the environment acts on the system through 
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classical stochastic fields. The effect on the system is the same as that of other environments, 
provided a suitable ensemble average is taken. The coupling to the system occurs through 
spin operators. We distinguish two systems, in one of which the coupling operator commutes 
with the system Hamiltonian, in the other it does not. For both systems, we implement a 



Grover search algorithm 23|, [24| and demonstrate the effect of different environments. 



II. SYSTEM AND ENVIRONMENT 



A. System Hamiltonian 



We use a system of two qubits to compare the effect of different colored noise for different 
energy level structures. The Hamiltonian of the system is 

1 



h s = -h 

2 



lli 2 2 22, j 1 2 

j z a z + u z a z - u x a x + irJa z a z 



(1) 



Here a % x , z denotes the x or z component of the Pauli matrix for spin i, uj x / z describes the 
strength of the magnetic field along the x or z axis, respectively, and J denotes the coupling 
constant. 

For suitable parameter sets, this Hamiltonian can execute a controlled NOT (CNOT) op- 
eration in a single step, without external control operations jl]]; for a wider set of parameters, 
two-qubit gate operations can be executed that fall into the CNOT equivalence class, i.e. 
they are equivalent to the CNOT operation up to single qubit operations. Table [J summa- 
rizes two parameter sets: System I ( left hand column ) gives the parameters for the CNOT 
operation, system II ( right hand column) the parameters for a CNOT-equivalent operation. 
We write the Hamiltonians of the 2 systems H\ and Hy. Their eigenvalues are also shown 
in Table [D The transition angular frequencies are u nm = (E n — E m )/K for each system, 
where n, m = 1, 2, 3, 4. The different energy level structures of the two systems provides the 
possibility to suppress decoherence induced by the coupling to the environment. 

B. Coupling to the Bath 

We now consider decoherence processes that are induced by a coupling to the environment 
that takes the form 

H(t) = H S + h<Ks(t)A. (2) 



TABLE I: Parameters for the two systems in units of it J, where J is the coupling constant. 





System I 


System II 




0.378 


0.378 




1 


1 




2.272 


1.136 


E X 


-1.32ft 


-0.961ft 


E 2 


-1.32ft 


-0.758ft 


E3 


0.948ft 


0.379ft 


E A 


1.70ft 


1.34ft 



For our purpose, the system operator A may be either a\, a 2 , or a\ + cy 2 . The bath term s(t) 
of the coupling operator may be either a classical random field or a quantum mechanical 
operator; for our purposes, it will be sufficient to consider it a time- dependent magnetic field 
with zero mean. 

For the random perturbation s(t), we consider stochastic functions with a Lorentzian 
spectral distribution of the power spectrum 

S(iU) = P + (a;-^ (3) 
and check the effect of the center frequency Uq on the decoherence rate in the two systems, 
where k describes the strength of S{ui). 

In the experiment, the coupling constant J had the value J = 215 s _1 . To simulate 
the effect of the environment, we generated the stochastic functions by digital filtering 
of a random time series. Each time series had a duration of 24.35 ms consisted of 80 
segments. For each time series we performed an experiment and summed over the individual 
experimental data. Figure [1] shows the power spectra of some time series s(t) and the RMS 



spectral density yS(u). Figure 2 shows the spectral density functions [25] for the four 
reservoirs that we compare, in relation to the transition frequencies of the system. 
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III. GROVER SEARCH IN NOISY SYSTEMS 



The elementary gates in the Grover search are the Walsh- Hadamard transform and the 
controlled phase reversal I\ x ) where \x) denotes a computational basis state, e.g. 1 00) or 
1 11). I\ x ) can be implemented by CNOT gates and one- qubit operations [19]. Hence we 
first implement the CNOT gate using the evolution under H s (pQ). 



A. Implementation of CNOT gates 



As stated in the introduction, the system Hamiltonian generates a CNOT-equivalent 
operation without additional gate operations, i.e., for a suitable time tc, U(t) = e ~ ltcHa ^ h 
becomes 

C e = CRl((j)), 



where R\{4>) = e 1 ^! 2 and 



C 



\ 



(4) 



10 
10 
1 
10 

We diagonalise the Hamiltonian H s = VDV^ . Here, V represents the eigenvector matrix 



V 



v 



( a x a 2 ^ 

Pi /3 2 

7i 72 

5i 5 2 



(5) 



and D is the diagonal form of the Hamiltonian, with eigenvalues 



K2 = ^[vl±Muif + {ui+ J -y] 



A 



3,4 



j , 



To determine the required evolution time tc*, we choose the target operator as C e and 



calculate the fidelity 



26] 



,t) = \Tr[U(t)Ct)\/4. 
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F((f),t) = \ e ~ itXl/h + e - itX2/h - ^ ( e -itA,/s _ e -itA4/A)j/ 4> ( 6 ) 



We find 



:^ 2 ) 2 + k - 1) 2 

Numerical solutions for J = 215 s^ 1 are for system I (jr c — 0, = 6.15 ms, resulting in 
F > 0.9999, and for system II <p]l = 0.18vr, t]l = 4.05 ms, or 0g = 0.52vr, = 12.18 ms, 
F > 0.999. Figure [3] shows the time dependence of the fidelity for the three cases. 



B. Grover search 

Using CNOT gates, one obtains the controlled phase reversal 1^ = W 2 CW 2 and /| o) = 
e i(n/2)a z j where W 2 denotes the Walsh-Hadamard transform for qubit 2. The single-qubit 
gate operations required for the Grover algorithm are implemented by short radio-frequency 
pulses, whose duration is negligible compared to the two-qubit gate. It is therefore sufficient 
to consider the coupling to the reservoir during the evolution under H\ and Hi 1 . 

Including the perturabations, the total Hamiltonians are H\{t) = H s + hiraSk{t)A and 
H\}{t) = Hi 1 + hnask(t)A, respectively, where Sk(t) has been normalized and a represents 
the strength of the perturbation. The perturbed Hamiltonian generates a perturbed phase 
reversal 1^ v x \, which deviates slightly from the ideal operation L x \, and which differs for each 
instance of Sk{t). 

The initial state for the Grover search is the uniform superposition = ( 1 00) + |01) + 
1 10) + |ll))/2 obtained by applying W 1 ' 2 to 1 00) . We choose the target state as |11). 
Figure H] shows the sequence of gate operations for the full Grover search algorithm G = 
W 1 ' 2 I\qq)W 1 ' 2 I\ii) where we have used W = W~ x . 

In system I the a single CNOT gate takes 6.09 ms, in the second system either 4.05 or 
12.18 ms. To make the duration of the algorithm in both systems comparable, we replaced 
the second CNOT operation in system I by CNOT 3 , which takes 18.26 ms. Since the 
algorithm includes 2 CNOT gates, the total duration is close to 24 ms in both systems. 

In each experiment, we start from the pseudo-pure state |^o)- The perturbed Grover 
search Gk = W 1,2 I\qq)W 1,2 I\ii) transforms it into \^k) = Gk\^o) and the corresponding 
density matrix into pk = l^^i^kl- Averaging over the individual signals gives the average 
(mixed) density matrix 
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28] 



1 M 

P=uT,Pk- (7) 

1V1 k=i 



C. Decoherence during the search process 



To quantify the loss of coherence by the environmental perturbation, we measured the 
purity of the Grover search process by averaging the purity of the final states for proper 
input states. For this purpose, we chose a set of states that is uniformly distributed over 
the Bloch sphere. The uniformly distributed set of input states consists of the 36 states 
l*S?> = \1>a)\1>b), K b = 1, 2, . . ., 6) where |V> a ,&> e {|0>, |1>, (|0> + (|0> - \l))/y/2, 

(|0)+i|l))/V2, (|0)-i|l»/>/2}. 

The average purity for the algorithm is then P = ^ Y^=i Tr{[p^} 2 } where p^ denotes 
the output density matrix after completion of the quantum search for the input states 
p^ = \^in^)(^in\- Note that |\I/o) is one of the 36 input states. 

When the systems are embedded in the reservoirs shown in Figure [2J the final states p^ 
can be calculated by solving the Bloch- Redfield equations. In the eigenbase of the system 
Hamiltonian H s in Eq. ([T|), the Bloch- Redfield equations are 



Pnm 'i^nmPnm ^ ' RnmklPkh 

k,l 



(8) 



where R nm ki denotes the partial decoherence rates 



R 



nmkl 



film J]] h-nrrk + $nk J]] -^Irrm ^linnk ^kriml- 



(9) 



Mmnk denotes the element of the relaxation tensor. Eq. Q shows that one can obtain R nm ki 
through the real parts of the relaxation tensor represented as 



Re{Mmnk} — —S(uj nk )(Ai m A nk ). 

47T 



(10) 



When A = a z , it takes the form 



A 



0-10 

-0.6606 0.7507 

-10 

0.7507 0.6606 



A 



n 



^ -0.8695 0.4940 ^ 
0-10 
-10 
0.4940 0.8695 



V 



;n) 



in the energy representation of the systems I and II, respectively. When A = al, it is 
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diagonal because [cr*, H s ] = 0. The matrices in the two systems are then 



A 1 



f -1 0^ 



10 
0-10 
1 



V 



.4 
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^10 0^ 
0-100 
0-10 
1 



(12) 



Eqs. (15UTU1) show that decoherence can be suppressed by choosing the parameters of H s 
such that the elements of the tensor R get small Qualitatively, the influence of the 
environment depends on the size of \S{ui n k)A n k\. 



D. Numerical evaluation 



We first present a numerical evaluation of Eq. (jSJ) to calculate the purity for the 
two systems. Figures [5] and [6] summarize the result. Figures [5] (a-c) show the purity for the 
perturbation operator A = a\ and the reservoirs R1-R3. In Figure (a), the increase of a, 
has very little effect on P 2 because S{u l ^ k ) « 0. However, P 1 decreases significantly, because 
S{uj\ 2 )A\ 2 = 0.7507. Comparing the two systems in Rl, one finds that system II is more 
robust. In Figure (J5] b) the situation is reversed and system I is more robust. In Figure (jSl 
c), Pi decreases faster than P 2 because ^(u;^)^^! > ^(a^-AjjjJ. These results illustrate 
the possibility to suppress decoherence by choosing an appropriate energy level structure. 

When A — &1, it is diagonal in the eigenbase of the Hamiltonian. In this case, the 
reservoir does not induce transitions, but only causes dephasing, according to Eq. ATS"]) . The 
energy level structure has then only a small effect on the decoherence rates. ^From Eq. ( TTTJT) . 
one finds that the reservoir affects the quantum system only through its static part S(0). 
Figure [6] shows the resulting purities Pi and P 2 for the case that qubit 1 is coupled to the 
reservoirs R3-R4. In Figure (a) both Pi and P 2 remain close to 1 because S(0) ~ 0. In 
Figure (b), however, S(0) = 1, and both systems are affected in a similar way. 
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IV. EXPERIMENTAL PROCEDURE AND RESULTS 



A. Implementation of Hamiltonians 

For the experimental implementation, we chose Carbon-13 labelled chloroform (CHCI3) 
dissolved in d6-acetone as the quantum register. We chose the carbon as qubit 1 and the 
proton as qubit 2. The noise term in Eq. ([2]) is introduced by an offset variation of the 
transmitter on channel 2. The Hamiltonian of the two qubit NMR system is thus 

111 

H NM R,k{t) = -xhu z al + -huj 2 z a 2 z + -Tin J a l z a 2 z + %its k {t)A, (13) 

where A = a 2 or a\ and the coupling is J = 215 Hz. 

The transverse field in Eq. ([1]) is applied as a radio frequency (rf) field, which can be 
written as 

Hrf = -\huy x (14) 

in the rotating reference frame. 

In the experimental implementation, the perturbation and thus the total Hamiltonian are 
piecewise constant for short periods r = 304.38 fis. For each of these periods, we realized 
the total evolution e~ lTHk ^^ n as e~ lTHNMR - k ^ n e~ lTHr ^ n , i.e. by a short free precession period 



followed by a small flip-angle pulse |30|]. This is a good approximation when r <C 2ir/uj nm . 

For the experimental implementation of the reservoirs R1-R4 shown in Figure El we used 
M = 12 noise traces. The corresponding spectral functions are similar to those represented 
in Figures [2] (a-d). 

'he experiments start with the effective pure state 1 00) (00 1 prepared by spatial averaging 
32j. The pulse sequence [a} 2 x -[grad} z -[n / A} 2 X - ^-\n]l' 2 - ^-[-n]l' 2 -[-n / A] 2 -[grad] z 



31 



transforms the system from the equilibrium state to the effective pure state 1 00} (00 1 . Here 
a = arccos(27i/72) ~ tt/3, where 71 and 72 denote the gyromagnetic ratios of 13 C and 1 H, 
respectively, and [grad] z denotes a gradient pulse along the z-axis. denotes the evolution 
caused by H NMR for a time jj. The pulses are applied from left to right. The complete 
pulse sequences for the implementation of the Grover search in the two systems are shown 
as Figures [7] (a) and (b), respectively. 
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B. Grover Search in systems not coupled to reservoirs 



When the system is in |00)(00|, we experimentally measured the density matrix shown as 
Figure [S] through state tomography 33j. In either system the Grover algorithm is repeated 
12 times, using different noise traces. The target state is |11). The final NMR signals are 
obtained by summing the 12 signals acquired via the readout pulse. 

We first implement the Grover search in the systems without engineered noise, i.e., the 
noise signal asuif) {k = 1, 2, • • •, 12) is not applied to the quantum systems. After the 
completion of the search algorithm, the density matrices of the two systems are shown as 
Figures M (a-b), respectively. 

In the experiments, the imperfections of the rf pulses and natural decoherence cause 
errors in the search results. In order to distinguish these errors from those that are due 
to the engineered "noise" that we investigate here, we use the results of in Figures [9] (a-b) 
as the references for subsequent experiments in the systems coupled to the reservoirs. We 
denote these reference states as p\ and p Q . 

To estimate the effects of the errors caused by imperfections of the rf pulses and natural 
decoherence, we compare the experimental results to simulation data, where the rf pulses 
are perfect and no natural decoherence exists. These results are shown as Figures [9] (c-d) 
corresponding to (a-b), respectively. The overlap between p an d its corresponding simulated 
result is 0.97 for system I and 0.91 for system II. 



C. Search results in systems coupled to reservoirs 

To simulate the noisy reservoir and observe the decoherence effect, we apply M — 12 
different noise traces ask{t) to the system during the implementation of the Grover search 
and add the resulting signals. We quantify the resulting decoherence by the fidelity F p = 
Tr{p p). 

Figure [10] shows the resulting density operators for the case where the coupling operator 
is A = a\ and the reservoir is R1-R3. The upper row corresponds to Fig. I, the lower row 
to Fig. II. In all cases, the coupling strength was set to a = 63.66 Hz. The fidelity F p 
shown in the figures was calculated as the overlap between the states resulting from the 
noisy experiment and those from the experiment without the reservoir. 
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The observed results are in good agreement with the predictions from the numerical 
simulations shown in Figure [51 For example, the reduction of P\ in Figure [5] (a) leads to the 
low fidelity F p = 0.69 in Figure [TU] (a), while F p in Figure [TD] (d) is 0.92 and P 2 in Figure [5] 
(a) remains close to 1. In Figures [TD] (c) and (f), the fidelity in system II is larger than that 
in system I, in good agreement with the result that P 2 > P\ in Figure (c). 

Figure [11] shows the search result for the case where the coupling operator is A = o\ 
and the systems are coupled to reservoirs R3-R4. In Figures (a) and (c), a = 63.66 Hz; in 
Figures (b) and (d), a = 25.46 Hz. The experimental results agree with the results in Figure 
El The much higher fidelity of the first column show clearly that the reservoir affects the 
system only through S(0). 



V. GENERALIZATIONS 



The above description of the decoherence process uses the semiclassical approximation, 
where the environment interacts with the system through classical fields. The results are 
easily generalized to the case of a quantum mechanical environment. For this purpose, we 
describe the total system (quantum register plus bath) by the Hamiltonian 

H tot = H S + H B + H I (15) 

where Hb denotes the Hamiltonian of the bath, and Hi denotes the coupling between the 
system and the bath. For the purpose of comparison we choose Hj = AX where X denotes 
an operator of the bath. 

In the quantum mechanical description, the dynamics of the quantum register are ob- 
tained by tracing over the degrees of freedom of the environment. It is thus possible to 
recover the Bloch- Redfield equations. In the eigenbase of H s , the rates A/ mn fc flHl fTOj) become 



lalisl 



Aimnk = A lm A nk / -^e- luJ " kt (X(t)X(0))dt (16) 
Jo n 

where X(t) = e %H Bt/%Xe~ iHBt ' h . The brackets (. . .) denote the thermal average over the bath 
degrees of freedom. The Fourier transform of this correlation function corresponds to the 
spectral function S{oj n k) in Eq. (jTOl) . This means that the fully quantum model described 
by Eq. f[T5|) can be mapped to the quantum system under classical noise described by Eq. 
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(J2J) [28]. Consequently our results are equally applicable to the fully quantum- mechanical 
case. 

Besides the coupling between the quantum system and its environment, pulse imperfec- 
tions (i.e. nonideal gate operations) also induce decoherence. Our results can be easily 
generalized to investigate the effect of the pulse imperfections on the decoherence rate. The 
generalization is illustrated by rewriting Eq. (J2]) as 

H(t) = \h [nl(t)al + u\al - u£o* + vJo\e^ (17) 

where fll(t) — uj\ + 2irs(t), when A — a\. f^(i) denotes the strength of the pulse that 
randomly fluctuates about a;*, and 2iTs(t) describes the fluctuation. Using our methods, one 
can discuss the effects of different fluctuations and search for experimental conditions that 
minimize the effect of pulse imperfections. 

VI. CONCLUSION 

We have investigated, experimentally and theoretically, the effect of different reservoirs on 
the decoherence of quantum registers during the execution of a quantum algorithm. While 
we have used a semiclassical system for these investigations, the results are easily adapted 
to a quantum mechanical environment, such as a spin-boson model. 

The system Hamiltonian as well as the coupling operator determine whether the envi- 
ronment causes pure dephasing or also induces transitions. The situation that is probably 
most relevant for quantum information processing is the case where the coupling operator is 
diagonal in the eigenbase of the system Hamiltonian. In this case, the environment causes 
pure dephasing and only the static part of the perturbation, oc 15(0) | causes decoherence. 

While we have chosen a 2-qubit system for this investigation, the results are completely 
general and can be applied directly to multi-qubit systems. It is possible to use this method of 
simulating dissipative quantum systems for related phenomena, such as dissipative quantum 



phase transitions 



341 ] . In the field of quantum information processing, our results indicate 



possible ways for suppressing decoherence. 
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FIG. 2: (Color online) Normalized spectral functions for the four reservoirs (R1-R4) with T = 100 
s _1 and uio = 0%, uj^, (0% + ujI[)/2 w wlf and 0, respectively. The thick curves represent 

the generated spectral functions for each reservoir, which were generated by averaging over 2500 
noise signals. The thin curves represent the corresponding theoretical spectral functions. The 
transition angular frequencies of the systems I and II are marked by "*" and "+". 
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FIG. 3: (Color online) Time dependence of the fidelity for the CNOT-equivalent operation. In 
system I, shown as the solid curve, the proper evolution time is 6.15 ms, and 4>c = 0. In system 
II, shown as the dash-dotted and dashed curves, the proper evolution times are 4.05 ms and 12.18 
ms for 4>c = 0.18-7T and 0.527T, respectively. 
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FIG. 4: (Color online) Gate sequence for the Grover search in an open quantum system. The time 
order is from left to right. W denotes the Walsh-Hadamard transform, and the red arrows denote 
the interaction between the system and reservoir. The duration of the gates W and e lna ~z/ 2 is so 
short that it can be ignored. 
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FIG. 5: (Color online) Purity of Grover search process as a function of the coupling strength to 
the environment when qubit 2 is coupled to the reservoirs in R1-R3, respectively. The data points 
obtained in systems I and II are marked by "*" and " x", respectively. 
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FIG. 6: (Color online) Same as Fig. but for qubit 1 coupled to reservoirs R3-R4. 



19 



[-*] x [*/2] y UdyH, 



[-71/2] , 



[-7T/2] y H x [_ n/2 ] y [0.4871] x l-n/2] y [-0.52tt] „ 



C1 
|0> 

H2 
|0> 



C1 
|0>. 



[^2] y M x [e,] x [-*], 



[e„L [-^2] M x [e 2 ] x 



H2 
|0> 



20 Loops 



60 Loops 



40 Loops 



40 Loops 



FIG. 7: Pulse sequences to implement the Grover search in systems I shown in (a) and II shown in 
(b). The flip angles Q\ and 02 for the two systems differ by a factor of 2. During the delay denoted 
by r, the systems evolute under the natural Hamiltonian and the noise signal of Eq. f)13[) . 
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FIG. 8: (Color online) Experimentally measured density matrix when the system lies in the initial 
pseudo-pure state |00)(00|. 
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FIG. 9: (Color online) Experimentally measured density matrices after the completion of the 
Grover search in the systems I [shown as (a)], and II [shown as (b)] when the reservoir is not 
applied. The target state is chosen as |11). The matrices have been normalized. Only the real 
parts of the elements are plotted. The imaginary parts are less than 18%. In order to estimate 
the errors caused by the imperfections of rf pulses and natural decoherence, Figures (c-d) show 
the simulated results obtained by NMR simulator where the rf pulses are perfect and no natural 
decoherence exists, corresponding to Figures (a-b), respectively. 
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FIG. 10: (Color online) Experimentally measured density matrices after the completion of the 
Grover search when qubit 2 is coupled to reservoirs R1-R3, shown as the three columns from left 
to right. The two rows of figures show the results obtained in systems I and II, respectively. F p 
denotes the fidelity of the search result with respect to the corresponding noiseless result. 
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FIG. 11: (Color online) Experimentally measured density matrices after the completion of the 
Grover search when qubit 1 is coupled to reservoirs R3-R4, shown as the two columns from left to 
right. 
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